2009 Fermi Symposium, Washington, D.C., Nov. 2-5 

Cross Calibration of Imaging Air Cherenkov Telescopes with Fermi 

Manuel Meyer, Hannes-Sebastian Zechlin, Dieter Horns 

Institut fiir Experimentalphysik, University of Hamburg 
Lumper Chaussee 149 D-22767, Hamburg, Germany 
(Dated: August 27, 2010) 

An updated model for the synchrotron and inverse Compton emission from a population of high energy 
electrons of the Crab nebula is used to reproduce the measured spectral energy distribution from radio to high 
energy y-rays. By comparing the predicted inverse Compton component with recent Fermi measurements of 
the nebula's emission, it is possible to determine the average magnetic field in the nebula and to derive the 
underlying electron energy distribution. The model calculation can then be used to cross calibrate the Fermi 
observations with ground based air shower measurements. The resulting energy calibration factors are derived 
and can be used for combining broad energy measurements taken with Fermi in conjunction with ground based 
measurements. 



L INTRODUCTION 

The Crab Nebula is probably the best studied object in 
astrophysics (for a recent review see e.g. Hester lfT6l ). It 
is the remnant of a core-collapse supernova which occured 
in 1054 AD at a distance of d w 2kpc 11271 - Near its ge- 
ometric center resides a pulsar which continuously injects 
a wind of ultra-relativistic particles into the nebula. The 
wind terminates at a shock front where electrons are pitch 
angle isotropized, forming a broad power-law in energy. 

Observations of the synchrotron and inverse Compton 
nebula have been carried out in every accessible wave- 
length resulting in a remarkably well known spectral en- 
ergy distribution (SED). The large area telescope (LAT) 
onboard the Fermi satellite [26] has recently measured 
the flux of the y-rays of the Crab between 100 Me V and 
300 GeV with unpreceded accuracy. 

These observations can be used for a cross calibration 
between the Fermi/LAT and ground based air shower ex- 
periments. This leads to the elimination of the systematic 
uncertainties on the absolute energy scale of typically 15 % 
for imaging air Cherenkov telescopes (IACTs). 

For the cross calibration we use our approach presented 
in Meyer et al. E4l : based on the work of Hillas et al. 
Ifl8l , we extract the distribution of electrons from the 
synchrotron spectrum under the assumption of a constant 
magnetic field strength throughout the nebula. Using this 
electron distribution in conjunction with the seed photon 
fields as extracted from observations, we obtain a detailed 
prediction for the inverse Compton emission. 

The average magnetic field inside the nebula is deter- 
mined by fixing the model to the Fermi/LAT observations. 
The model can directly be used as the basis for the cross 
calibration. For this purpose energy scaling factors are de- 
rived to correct the measurements of the ground based in- 
struments to the model. 

The article is organized as follows. In section [TTJ the 
model for the SED is quickly reviewed, where the details 
and the discussion of the underlying electron spectrum are 
spared out and we refer the reader to Meyer et al. [24 1. 
The results of the cross calibration are presented in section 
llirjtogether with an application to extract limits on the dif- 



fuse y-ray background at TeV energies. Section IV sum- 
marizes the article and gives a short outlook on possible 
future work. 



II. MODEL OF THE SPECTRAL ENERGY 
DISTRIBUTION OF THE CRAB NEBULA 

A population of relativistic electrons is assumed to be 
distributed in a spherical volume following a Gaussian 
density distribution with its maximum at the nebula's cen- 
ter. The width of the Gaussian is parameterized in order 
to reproduce the shrinking size of the nebula with increas- 
ing frequency. The volume of the nebula is assumed to 
be filled with an entangled magnetic field of constant field 
strength. Within the volume occupied by the electrons, 
various seed photons are upscattered. The effective den- 
sity of seed photons « see d is simply found by convolving 
the electron density with the photon density [see Q~8| for 
further details]. The total emission from the nebula is then 
found by integration: 



J d r n(y)(X^ + 4 c ), 



(1) 



with X.v the single particle emission functions for syn- 
chrotron (Sy) and inverse Compton processes (IC) [see e.g. 
0: 
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where we have averaged the pitch angle to give y/2/3 and 
cr T denotes the Thomson cross-section. The critical fre- 
quency v c is defined as 

3e 
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and K^/j(x) stands for the modified Bessel function of frac 
tional order 5/3. Introducing the kinematic variable q, 

hv 
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(5) 



eConf C091122 



2 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 



the IC distribution function fic can be written as 



fic(e,v,y) 
2q\nq 



+ (l+2q)(l-q) 



1 [4ey^/(mc 2 )] 

2 1 + Aeyq 2 l(mc 2 ) 



d-q)- (6) 



For the inverse Compton channel photons from several 
seed photon fields are taken into account: (1) synchrotron 
radiation, (2) emission from thermal dust, (3) the cosmic 
microwave background (CMB), and (4) optical line emis- 
sion from the nebula's filaments. The seed photon density 
n see d in Eqn. [3] is the sum of all these components. Like 
the electron population, the spatial photon densities are ap- 
proximated by Gaussian distributions whereas the photon 
density of the CMB is taken to be constant throughout the 
nebula. The spatial variance of these distributions is also 
energy dependent. The resulting broadband energy distri- 
bution is shown in Fig. [T] The electron spectrum is varied 
until the resulting synchrotron spectrum reproduces the ob- 
servational data. 

The compilation of data used here is summarized in Aha- 
ronian et al. 12 and references therein. Additionally, new 
data are added which are listed in Table Q] The solid black 
curve in Fig. |T| is the sum of all contributions including 
synchrotron and IC emission as well as thermal emission 
from dust in the nebula and optical line emission from 
the filaments. For the thermal dust emission a graybody 
spectrum is used and a temperature of T — 93 K was de- 
rived by fitting the combined spectrum (thermal and non- 
thermal emission) to the data (solid gray line in Fig. [TJ. 
The flux of the line emissions (orange solid line) is taken 
from Davidson ifTTIl . Davidson & Fesen lfT2l and Hester 
et al. IfTTIl . The optical line emission of the filaments in 
the nebula is estimated in the following way: The high res- 
olution spectral observations of individual filaments have 
been corrected for extinction [ 17 1 and scaled to match the 
global emission from the filaments [see e.g. the discussion 

inim. 

The FIR observations from Spitzer, ISO, and Scuba (or- 
ange and magenta circles in Fig. [TJespectively) deviate 
from the power-law extrapolation of the radio spectra. In 
the framework of two distinct electron populations, the 
shape of the continuum is naturally explained by the tran- 
sition of the two synchrotron emission components. 
The dashed blue curve indicates the total synchrotron and 
IC emission which results from the contributions from 
wind electrons whereas the dashed red line shows the con- 
tribution of the radio electrons. The dashed cyan lines 
show the total synchrotron and IC flux from both popu- 
lations, respectively. 

The electron spectrum was adapted such that the result- 
ing synchrotron emission reproduces the power law mea- 
sured with XMM-Newton. Both, the XMM-Newton and 
INTEGRAL (with the instruments SPI and IBIS/ISGRI) 
observatories are calibrated on the basis of detailed simu- 
lation and laboratory measurements. This approach differs 
from widely used corrections of the response function in 



Energy Band Instrument Reference 

Sub milimeter ISO & SCUBA Green et al. (TJ 
to far infrared SPITZER Temim et al. 1251 
X-Ray to XMM-Newton Kirsch et al. (23) 
y-rays SPI Jourdain & Roques 1 19 1 

IBIS /ISGRI Jourdain et al. |20l 
Fermi / LAT The Fermi collaboration 1 26 1 
VHE H.E.S.S. Aharonian et al. (3) 

MAGIC Albert et al. 

TABLE I: References for the observations used for the SED. All 
other data are taken from Aharonian et al. |2| and references 
therein. 



order to re-produce a specific spectral shape and flux of the 
Crab nebula. The corrected measurements are therefore 
model-dependent and have not been included here. Fur- 
thermore, XMM-Newton is able to spatially resolve the 
Crab, whereas the measurement by SPI and IBIS/ISGRI 
may include contributions from the pulsar, possibly lead- 
ing to higher fluxes in comparison to the XMM-Newton 
observations. Note, that the difference in flux normaliza- 
tion between XMM-Newton and SPI are beyond the sys- 
tematic errors quoted. The model calculations shown here 
are fixed to the XMM-Newton spectra which then natu- 
rally underpredict the SPI measurements. The shape of 
the spectrum measured by SPI has been taken into account 
as both the power-law below 100 keV as well as above 
smoothly connect to the spectra measured at the low en- 
ergy end by XMM-Newton and with Comptel at higher 
energies. 

The predicted IC emission and the Fermi observations 
are used to determine the average magnetic field. This is 
shown in Fig. [2] The IC fluxes due to the different pho- 
ton fields and electron populations add up to give the total 
black solid line. A standard x 2 -minimization is used to 
determine the best value for the average B-field. Since a 
varying magnetic field also changes the synchrotron flux, 
the underlying electron spectrum is varied accordingly to 
compensate for the change, i.e. the synchrotron flux re- 
mains constant. Taking the systematic energy uncertainty 
on the global energy scale of the Fermi data into account, 
AE/E = (sys.) [see e.g. [T), the average B-field is 
found to be 



B = (l24 + 6 (stat.) (sys.)) //G. 



(7) 



The discussion of this result in the light of MHD calcu- 
lations ETI can be found in Meyer et al. 11241 . It is worth- 
while noting that the magnetic field derived here is less 
than half the value of the commonly used 300 jiG. 



A. The electron spectrum of the nebula 

The underlying electron spectrum dN e j/dy is the crucial 
quantity that determines the shape of the SED [see 24., for a 



eConf C091122 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 



3 



log 10 (E/eV) 

-5 5 10 




10 15 20 



log 10 (v/Hz) 

FIG. 1: Broadband SED of the Crab nebula. See section [TT] for details. The two black curves correspond to the models described in 
Meyer et al. [ 24 1 
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FIG. 2: The total IC flux due to different seed photon fields and the Fermi data points are shown. For the numbering see the text. 



detailed discussion]. It consists of the two afore mentioned 
electron populations where the radio electron spectrum is 
given by 

dN r [Ky- Sr for y[<y<y> 2 , 

= i ■ (8) 

ay 

[ otherwise, 



The radio electrons were probably injected in the phase of 
rapid spin-down during the initial stages of the pulsar-wind 
evolution [6]. The values for yL y' v and y' 2 are summarized 
in Table [II] Above y r 2 and below y\ a sharp cut-off for the 
spectrum is chosen, i.e. dN e \/dy = Ofory < y\ and 7 > y r 2 . 

The wind electrons produce via synchrotron emission 
the bulk of the observed SED above sub-mm/FIR wave- 
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Population No lny In y [ lny 2 S 
Radio 120.0(1) - 3.1 12.1(7) 1.60(1) 
Wind 78.6(3) 19.5(1) 12.96(3) 22.51(3) 3.23(1) 

TABLE II: Cut-off and normalization energies together with the spectral indices S of the electron spectrum. See text and Eqn. [9] for 
further details. 



lengths. The wind electrons are constantly injected down- 
stream of the wind shock (hence the name). The radia- 
tively cooled spectrum of the wind electron has a spectral 
index of S w = 3.23 = 2.23 + 1 which can naturally be 
explained by ultra-relativistic 1 st order Fermi acceleration 
with synchrotron cooling [see e.g. [22). 
An additional feature present in the hard X-ray spectrum 
which follows a broken power-law with a break at ~ 70 kev 
requires a break with AS - 0.43 in the electron-spectrum. 
We tentatively relate this break to the injection mecha- 
nism, given that it can hardly be related to energy depen- 
dent escape (the X-ray emitting electrons suffer cooling 
well before escaping the nebula). The value of AS could 
hint at an energy dependent effect similar to diffusion in a 
Kolmogorov-type turbulence power spectrum. 
For high and low energies the wind electron spectrum cuts 
off super-exponentially, 
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where the values for y^ , y™ and y™ are listed in Table [TTJ 
The total energy of the radio and wind electrons, respec- 
tively is found to be 

X°° d^V' 
y-^-dy = 3.10 x 10 48 ergs, (10) 

E w = mc l f y — - dy = 2.28 x 10 4S ergs, (11) 
Ji dy 

indicating that the total energy is much smaller than the 
integrated energy released through the spin-down of the 
pulsar. The fact that both relic electrons and wind electrons 
share roughly equal energy is probably coincidental. 



III. CROSS CALIBRATION OF IACTS & FERMI 

The updated model of the SED of the Crab Nebula 
provides an opportunity for the cross calibration between 
ground based air shower experiments and the Fermi large 
area telescope. The method is demonstrated here with 
the imaging air Cherenkov telescopes (IACTs) HEGRA, 



H.E.S.S. and MAGIC but is applicable to any other ground 
based air shower experiment. 

The energy calibration of IACTs is done indirectly with 
the help of detailed simulations of air showers and the de- 
tector response. However, the remaining systematic uncer- 
tainty on the absolute energy scale of typically 15 % leads 
to substantial differences in the observed flux and position 
of cut-offs in the energy spectra between different IACTs 
and also between Fermi/LAT and IACTs. 
Since the observed energy spectra are usually quite broad 
in energy, the position of features in the spectra are not 
useful (and may be time-dependent for some objects) for 
cross calibration. On the other hand, cross calibration be- 
tween Fermi/LAT and IACTs provides indirectly a means 
of benefitting from the careful beam-line calibration of the 
Fermi/LAT [see e.g. 8|. For this reason, the average mag- 
netic field used in the model was fixed to Fermi observa- 
tions. 

The cross calibration is now accomplished in the follow- 
ing way: for each IACT an energy scaling factor stact is 
introduced such that 



E = E ■ iiACT- 



(12) 



The scaling factor siact is determined via a x 1 ' 
minimization in which the energy scale is changed accord- 
ing to the formula above until the data points reproduce 
the model best. The scaling factors for the different instru- 



ments are listed in Table III together with the statistical 
errors and the reduced x 1 values of the fit. The statisti- 
cal uncertainties were obtained by summing the errors of 
the^ 2 -fit and the statistical errors of the model in quadra- 
ture. The latter are mainly due to the uncertainties on the 
B-field of Eqn. [7] To illustrate the result, Figure |3(a)| and 
[3(b)] compare the unsealed data points of the Crab nebula 
with the scaled ones. It is evident that the data points fit the 
model better after scaling. All scaling factors lie within the 
afore mentioned 15 % energy uncertainty of the IACTs. 

The application of the cross calibration eliminates the 
systematic uncertainty of the energy scale of the IACTs 
and adjusts the energy scale to the one of the Fermi/LAT. 
Since the model relies on the Fermi/LAT measurements, 
the Fermi/LAT's absolute energy uncertainty remains. 
This implies an improvement for the systematic uncer- 
tainty on the energy scale of IACT measurements from 
±15 % to However, the model itself also con- 

-10 % 

tributes to the systematic uncertainties. The main uncer- 
tainty stems from the fact that a constant magnetic field is 
assumed which is very likely not the case. It should be also 
noted, that the cross calibration factors hinge mainly on 
the high statistics measurement of the IACTs at low TeV- 
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energies which are produced by electrons that co-exist in 
the same volume in the nebula as the electrons emitting at 
lower energies (i.e. in the Fermi/LAT energy range). 

As a first application of the cross calibration, we derive 
upper limits on the diffuse y-ray background. 

Both Fermi [1J, and H.E.S.S. l4l [T5l have measured the 
cosmic ray e~ + e + spectrum. Unlike Fermi, the tele- 
scopes from H.E.S.S. cannot accurately distinguish be- 
tween showers induced by electrons (or positrons) or pho- 
tons, such that up to as 50 % of the observed electromag- 
netic air showers could be induced by photons. Hence, 
H.E.S.S. actually measures electrons and diffuse back- 
ground photons. Taking the difference of the two measure- 
ments we can derive an upper limit on the flux of the y-ray 
background. The scaling factors derived above are now 
used to convert the IACT data to the same energy scale of 
Fermi, which reduces substantially the systematic uncer- 
tainty on the observed flux given that the electron spectrum 
follows a soft power-law with E~ 3 . 

The upper limits were derived by subtracting the two 
fluxes from the overlapping region of the measurements. 
This corresponds to the first six H.E.S.S. points of the 
low energy analysis in Figure |4] The remaining system- 
atic energy uncertainties, denoted by the green and yellow 
bowties, were taken into account for the derivation of the 
upper limits: the flux points of the H.E.S.S. measurements 
were shifted to their maximum value allowed by the sys- 
tematic uncertainties while the Fermi points were shifted 
to the minimum value. Hence, the result represents a con- 
servative approximation of the upper limits. 

An important result of the cross calibration is that the 
peak in the spectrum observed by ATIC [ 10 1 appears more 
unlikely after applying the scaling factors. 

IV. SUMMARY AND OUTLOOK 

An updated model for the SED of the Crab nebula has 
been introduced. It incoiporates a new electron spectrum 



that consists of two electron populations that have been 
studied extensively in the past [see e.g.|6l|7). 

The average B-field was derived by fixing the IC flux to 
Fermi/LAT measurements. 

The model makes it possible to derive energy scaling 
factors that eliminate the systematic energy uncertainties 
of ground based air shower experiments. An application of 
the cross calibration to the diffuse y-background has been 
presented and upper limits have been derived. Moreover, 
the excess measured by the ATIC collaboration seems un- 
likely with the scaled H.E.S.S. observations. 

Nevertheless, an improved model for the emission of 
the Crab is conceivable. Such a model could comprise a 
spatially varying B-field and hence the fact that electrons 
emitting different energies are exposed to different field 
strenghts. 

To conclude, a way to establish the Crab as a true stan- 
dard candle in y-ray astronomy has been presented and 
versatile applications to other observations of bright steady 
or pulsed sources (e.g. the Crab Pulsar, the Galactic Center 
other Pulsar Wind Nebulae, etc.) are imaginable. More- 
over, an application of the cross calibration will help to 
improve dark matter searches and constraints on the extra 
galactic background light. 
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Fermi - sys. error 
ATIC >- 

H.E.S.S - l.e.a. (scaled) y 
H.E.S.S. (scaled) h 
upper limit on diffuse photons 
PAMELA (preliminary) t 



FIG. 4: e~ + e + Spectrum reported by H.E.S.S. and Fermi. The 
cross calibration is applied, hence the uncertainty on the global 
energy scale is eliminated. 
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